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Motivated by the state of the art method for fabricating high density periodic nanoscale de- 
fects in graphene, the structural, mechanical and electronic properties of defect-patterned graphene 
nanomeshes including diverse morphologies of adatoms and holes are investigated by means of first- 
principles calculations within density functional theory. It is found that various patterns of adatom 
groups yield metallic or semimetaiiic, even semiconducting behavior and specific patterns can be 
in a magnetic state. Even though the patterns of single adatoms dramatically alter the electronic 
structure of graphene, adatom groups of specific symmetry can maintain the Dirac fermion behav- 
ior. Nanoholes forming nanomesh are also investigated. Depending on the interplay between the 
repeat periodicity and the geometry of the hole, the nanomesh can be in different states ranging 
from metallic to semiconducting including semimetaiiic state with the bands crossing linearly at 
the Fermi level. We showed that forming periodically repeating superstructures in graphene matrix 
can develop a promising technique to engineer nanomaterials with desired electronic and magnetic 
properties. 

PACS numbers: 81.05.ue, 73.22.Pr, 63.22.Rc, 61.48. Gh 



I. INTRODUCTION 

Propagation of electron waves through the honeycomb 
lattice attributes exceptional features to graphene.^ Con- 
duction of electrons within one-atom-thick layer with 
minute scattering makes the observation of quantum ef- 
fects possible even at room temperature.'-''^' Experimen- 
tal investigations have reported the observation of half- 
integer quantum Hall effect for carriers in graphene and 
possible magneto-electronic device applic at ions .1^ Most of 
the unique properties of graphene are related to its mono- 
layer lattice structure, linearly crossed tt bands at Fermi 
level with electron-hole symmetry. Recently, we showed 
that the honeycomb structure with linear band crossings 
at Dirac points is also common to Si and Ge.'^ 

In an effort to make semimetaiiic graphene suitable 
for electronic applications, it has been functionalized to 
generate band gaps. It was theoretically shown that it 
is possible to induce band-gap opening produced by the 
adsorption of atomic hydrogen on graphene by choosing 
specific adsorption periodicity.^ It was also experimen- 
tally shown that B- and N-doped graphenes can be syn- 
thesized to exhibit p- and n-type semiconducting prop- 
erties that can be systematically tuned with the dopant 
concentration .El Effect of hydrogenation and the transi- 
tion metal atom adsorption on the transport properties 
of graphene was also investigated theoretically.'^^ Using 
symmetry arguments and tight binding calculations it 
was shown that the periodic structure of defects (such as 
B and N impurities) on graphene can exhibit semimetaiiic 
and semiconductor behavior.l^Moreover, weak perturba- 
tion potential forming a large hexagonal lattice in a two 
dimensional electron gas was shown to lead a massless 
Dirac fermion Hamiltonian with linearly crossing bands 
at Dirac points. ^^1^^^ 



A majority of the current studies on graphene is 
devoted to its chemical modification to create deriva- 
tives with different structures and properties. So far 
three known derivatives of graphene have been suc- 
cessfully achieved in chemical reactions: graphene ox- 
ide (GO),'i^ graphane (CH^pHU and recently fluoro- 
graphene (CF).^'^^^ Although GO is a wide band gap 
material that is important for device applications, its 
atomic structure, wherein the carbon atoms are deco- 
rated with epoxides, alcohols and carboxylic acid groups, 
is not suitable for nanoscale manipulations. CH obtained 
by exposing carbon honeycomb structure to hydrogen 
plasma, is another example of the graphene-based chem- 
ical derivative. Upon the hydrogenation, semimetaiiic 
graphene is converted into an insulator. CF, the two di- 
mensional counterpart of teflon, is the most recent focus 
of graphene research. 

Much recently, the fabrication of large graphene sheets 
having high-density array of nanoscale holes, called 
graphene nanomeshes (GNMs),'22l has been the land- 
mark in controlling the electronic properties at nanoscale. 
Additionally, the formation of one-dimensional periodic 
Stone- Wales type defects producing metallic nanowires 
on graphene matrix has also been reported. These 
recent advances have made mesh configuration a con- 
trollable parameters to monitor physical properties of 
nanostructures.'^^ Earlier, interesting effects of peri- 
odically repeating holes in the electronic and mechanical 
properties of graphene nanoribbons were predicted from 
the first-principles calculations.l^Sl 

In this paper, we apply supercell method to reveal 
the electronic, magnetic and mechanical properties of 
graphene which is patterned by various adsorbates or 
holes. The atomic structure of all adsorbates and 
holes are obtained after extensive structure optimiza- 



tion. These periodically repeating superstructures or 
nanomeshes display properties which are rather differ- 
ent from those of graphene. We showed that not all pat- 
terns of adsorbates or hole with a 2D hexagonal lattice on 
graphene have linear band crossing, but only those which 
have a specific rotation symmetries. However, depending 
on the size of patterns or holes and the repeat period- 
icity, a GNM can be in different states ranging from to 
semiconducting including semimetallic with linear band 
crossing at the Fermi level. 



II. COMPUTATIONAL METHODOLOGY 

The present study revealed crucial effects of the point 
group symmetry of nanomesh on the resulting properties. 
Here, we start with a brief discussion of hexagonal sym- 
metry and apply simple tight binding model of 7r-orbitals 
to reve al the effect of lattice symmetry on the band 
crossing .1221211 Graphene has the space group P6/mmm 
and point group symmetry Dg/j. At the-F point, the 
group of the wave vector is isomorphic to the point group 
De/t- However, irreducible representation of the wave 
vector point group turns into T)2h and D^h at high sym- 
metry points M and K (or K ), respectively. It was 
shown that the tight binding Hamiltonian with nearest 
neighbor hopping parameter, t — 2.7 eV 
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well approximates the vr-bands of perfect 
graphene.!^^ Here (ci) is the creation (annihi- 
lation) operator of a tt electron at the lattice site i. The 
first term is the on-site energy of each carbon atom and 
equals to energy of the 2pz orbital. Energy eigenvalues 
of graphene and other two hypothetical crystal having 
square and hexagonal lattices with single atom in the cell 
are calculated and the contour plots their energy band 
gap in BZ are shown in Fig|T] For graphene, energy 
dispersion is linear at the vicinity of the K- symmetry 
(Dirac) points and the Fermi velocity, which is linearly 
dependent to nearest neighbor interaction |)arameter, 
can be given by the expression vp = 3td/2h.^^ 

Hypothetical square lattice has semimetallic band 
structure as in Fig. [ijb). In the reciprocal space the 
band gap is closed along the boundaries of squares. How- 
ever, differing from graphene, these bands do not have 
linear dispersion. The structure having hexagonal lattice 
with single 7r-orbital per unitcell is shown in Fig|ljc). 
Such a structure having a single 7r-orbital in each unitcell 
may also be realized by adatom saturation of one type 
(A- or B-type) carbon atoms of graphene. In this case, 
structure has three-fold rotation symmetry and hence six 
nearest neighbors. In this case, while saturation yields a 
dramatic change in the electronic structure, linear band 
crossings at Fermi level still occur at the points on a cir- 
cle of radius R = in BZ. As shown in Fig. U\c), 
this circle passes through the corners of hexagonal BZ of 
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FIG. 1: Contour plots for band gap. (a) Perfect graphene. (b) 
A crystal having square lattice with single atom in the unit 
cell, (c) A crystal having hexagonal lattice with a single atom 
in the unit cell. Band crossing occurs along the yellow/light 
contours on which the band gap becomes zero. Band gap 
takes its maximum value at brown/dark contours. 



graphene. The Fermi velocity is calculated to be in the 
order of 10^ m/s in the vicinity of band crossing points. 

While the tight binding model allowed us to under- 
stand the general features of the band gap in different 
lattices, it fails to account for the reconstruction and re- 
bonding near the defect. In the rest of the paper we per- 
form calculations from the first-principles to investigate 
various types of defects. To this end we carried out spin- 
polarized plane wave calculations'^^ '^^ using local density 
approximation (LDA)'^^ and projector augmented wave 
(PAW)'^^ potential. Patterns of defects are treated using 
supercell geometry, where a minimum of 10 A vacuum 
spacing is kept between the adjacent graphene layers. Ki- 
netic energy cutoff, Brillouin zone (BZ) sampling are de- 
termined after extensive convergence analysis. For the 
plane-wave basis set, the kinetic energy cutoff is taken 
to be /i^|k -I- Gp/2m = 500 eV. For partial occupan- 
cies Methfessel-Paxton smearing method^ is used. The 
convergence criterion of self-consistent field calculations 
is 10^*^ eV for total energy values. By using the conju- 
gate gradient method, all atomic positions and the size 
of unitcell were optimized until the atomic forces were 
less than 0.05 eV/A. Pressures on the lattice unit cell are 
decreased to values less than 1 kB. 



III. ADATOM PATTERNED GRAPHENE 
NANOMESHES 

In this section we show how a periodic decoration of 
graphene by adatoms modifies the electronic structure. 
Here we considered H, F, O and Mn adatoms, which are 
adsorbed at different sites and form a (4x4) supercells 
on graphene host matrix as shown in Figj2|^a)-(d). We 
did not consider the interstitial or substitutional decora- 
tion, since experimental studies treating various foreign 
atoms, such as Ijg 0^21, liHl, m, cP, CcP, Fe and 
G(i^, Au and Pt!^, revealed that these atoms prefer to 
be adsorbed at various sites on the surface of graphene. 
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FIG. 2: Spin-polarized energy band structure of a periodic patterns consisting of the (4x4) supercells each having a single 
(a) hydrogen; (b) fluorine; (c) oxygen; (d) manganese adatom. (e) A similar hydrogen pattern forming the (2x2) supercell on 
graphene host matrix allowing significant coupling between adatoms. (f-h) Periodic patterns of two, eight and six hydrogen 
atoms in the (4x4) graphene supercell, respectively. Isosurface of charge density of bands crossing near the ii'-point are shown 
by insets. For the sake of comparison, linearly crossing tt and 7r*-bands of perfect graphene host matrix are also superimposed 
in the band structures. The zero of energy is set at the Fermi level Ef- The band gaps are shaded in yellow. All the bands 
presented in BZ corresponding to the (4x4) supercell. 



but none of them is adsorbed at the interstitial sites, nor 
is substituted for carbon atom. 

Owing to the unpaired electron, a single hydrogen 
adatom adsorbed to the (4x4) supercell has a spin polar- 
ized, semiconducting ground state with a net magnetic 
moment of /Lt=l /x^. Upon the adsorption of a hydrogen 
atom band structure of graphene changes dramatically. 
Instead of linear crossing of tt and tt* -bands at Fermi 
level, dispersionless impurity bands occurs with 0.1 eV 
indirect band gap. Similar to hydrogen atom, the most 
favorable adsorption site for a F atom on graphene is 
the top site of carbon atoms. Upon the adsorption of a 
F atom, sp^ bondings of three C-C bonds below F atom 
are dehybridized and form tetrahedrally coordinated four 
sp^ type bondings. Differing from the decoration of H 
adatom, ground state is nonmagnetic and as a result of 
odd number of electrons F decorated graphene becomes 
metalized as shown in Fig. |2jb). The band crossing at 
the K-points does not occur, since F adsorbed at the top 
of C atom changes the sixfold rotation symmetry to the 
threefold rotation symmetry. 

Oxygen atom favors the bridge site between two un- 
derlying C atoms. Upon the adsorption of an oxygen at 
the bridge site two underlying C atoms become buckled 
by 0.36 A. C-C and C-0 bonds are calculated to be 1.51 
and 1.44 A, respectively. Resulting structure is a non- 
magnetic semiconductor with 0.63 eV direct band gap 
(Fig. [2|^c)). Valence and conduction band edges occurs 
between the K- and F-points. The adatom at the bridge 
site breaks the six fold rotation symmetry and hence hin- 
ders the linear band crossing. 

The situation is different in the case of Mn, which is 



adsorbed above the center of a hexagon in graphene ma- 
trix and induces negligible deformation. Only six nearest 
C atoms raise slightly to higher (0.02 A) position rela- 
tive to the plane. Localized, non-bonding Mn-3(i orbitals 
form flat bands near the Fermi level. On the other hand, 
the sixfold rotation symmetry is maintained even after 
Mn atom adsorbed at the hollow site above the center of 
hexagon. Accordingly, the metallic structure with a net 
magnetic moment of 3 /is per cell allows linear crossing 
of graphene bands at the i^-points in Fig. |2jd). 

Having discussed the effect of periodic decoration by 
single adatoms, we next consider the periodic patterns 
of adatom groups. In Fig. [2]^e) we show electronic 
band structure corresponding to relatively denser hydro- 
gen coverage (C:H=8). Such a nanomesh created by one- 
sided decoration of four hydrogen atoms in the (4x4) su- 
percell gives rise to a relatively dispersive bands and a 
net magnetic moment of 4 /is per supercell. Since the 
six fold rotation symmetry of graphene is broken by ad- 
sorbed H atoms, linear crossing of bands do not occur. In 
a decoration involving two sides of graphene, where two 
adjacent C atoms of different sublattices are saturated 
from different sides as shown in Fig. |2](f). Since an equal 
number of A- and B-sublattice atoms are saturated, the 
structure is nonmagnetic semiconductor with a band gap 
of 0.8 eV. Another pattern derived from a graphane like 
domain consisting of 8 H atoms in Fig. |2jg) results in a 
band gap of 0.7 eV at F-point, but larger gap of 1 eV at 
the X-point. This nanomesh presents an electronic struc- 
ture rather different both from graphene and graphane. 
The electronic structure is, however, different for a pat- 
tern of six H atoms, which saturate six carbon atoms at 
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FIG. 3: Band structures showing the effects of the couphng 
between various patterns as a function of their size and the 
size of mesh (supercell). The linearly crossing of tt- and tt*- 
bands of graphene are shown by red/dashed lines. 



the corner a hexagon alternatingly from different sites; 
namely three of them adsorbed to A-sublattice from one 
side, remaining three adsorbed to B-sublattice from the 
other side. Even though sixfold rotation symmetry has 
changed to symmetry, both point group symmetries 
allow linear band crossing as seen in Fig. |2jli). This 
case demonstrates the crucial role played by the intrinsic 
symmetry of the pattern in determining the electronic 
structure .liilti^I 

Let us now take a closer look at the triangular pat- 
terns of adatom groups that have hexagonal symmetry. 
In Fig. |3] we plot these tt- and 7r*-bands in the vicinity 
(Ak=0.03 A~^) of if-symmetry point for different pat- 
terns on supercells of different sizes; namely the (4x4) 
and (8x8) (C128 and C32) supercells. Hg and H12 pat- 
terns which have 0.06 and 0.09 eV band gap at if-point, 
indicates that the band gap opening increases with in- 
creasing pattern size and hence increasing coupling. In 
the (8x8) supercell the interaction between periodically 
repeating patterns is hindered and hence linear cross- 
ing of TT- and 7r*-bands at ii'-symmetry point similar to 
bare graphene is attained. By expressing these bands as 
k = K -|- q and by neglecting the second order terms with 
respect to , the dispersion of the energy bands can be 
given as £'(q) ~ VFfi.|q|. Here the Fermi velocity vp 
calculated as 0.7x10® m/s (i.e. ^ 0.6 of the Fermi velocity 
calculated for Dirac Fermions in perfect graphene). 

We also present an example of Og pattern analogous to 
H12 decoration of (4x4) supercell (C32H12) of graphene 
in Fig. |3] In this pattern three O atoms are bound al- 
ternatingly to bridge site at one site, the remaining three 
to other site. Here it is seen that denser O patterns are 
not favored due to strong 0-0 repulsion. Linear crossing 
of bands at if-point can also be achieved by O adatoms 
forming a periodically repeating Og pattern. Small band 
gap in Fig. [3] can be closed if the supercell size is in- 
creased to hinder coupling between them. This example 
implies that the patterns similar to that done by using 
H atoms can be created by the passivation of p-orbital 
electrons with O atoms. In the case of fluorination of 
(4x4) graphene by Fg decoration, similar to Hg, linear tt 



and TT* -bands gets closer with very small (0.01 eV) band 
gap at if -point. Finally, the isosurface charge densities 
of these linearly crossing bands near if-point indicates 
that they mainly originate from graphene 7r-orbitals with 
small mixing from the adatom (see Fig. |2](h)). 



IV. HOLE PATTERNED GRAPHENE 
NANOMESHES 

Similar to adatom patterns on graphene, nanomeshes 
generated from the holes periodically patterned on 
graphene matrix exhibit also interesting features. This 
conclusion drawn from theoretical calculations are in 
line with the findings obtained from the fabrication of 
GNMs by means of block-copolymer assisted nanopat- 
terning process. It has been shown that GNMs hav- 
ing high-density periodic array of holes display promis- 
ing advantages relative to existing graphene devices .'^£1221 
Thanks to the advances in the p reparation of high qual- 
ity nanoscale hard maskJ^SHZMlin laboratory conditions, 
theoretical studies in this field become more relevant for 
applications. Here we carry out calculations for the holes 
having 1-2 nm repeat period and 2-10 A diameter. In 
Fig. |4]j^a) we describe the geometric parameters of C„ 
hole defects forming a hexagonal lattice, where n de- 
notes number of C atoms removed from graphene ma- 
trix to make a hole. For nanoholes, we define the hole 
size as the maximum diameter of the C„ hole defected 
region. After the creation of a Ci defect (i.e. single 
C vacancy), Jahn- Teller type distortion changes the po- 
sitions of surrounding C atoms slightly. The resulting 
structure attains a net magnetic moment of 1 [Ib- The 
origin of magnetism in defected graphene sheets and the 
character of electronic states induced by the vacancy re- 
sulting in fiat bands around Fermi level have been investi- 
gated by some recent studies .l^^^t^ Upon the removal of C 
atoms (n>l) graphene is reconstructed to result in a sig- 
nificant modification in atomic configuration around the 
hole. For example, after the relaxation of atomic struc- 
ture, C2 defected region becomes an octagon-shaped hole 
surrounded by 2 pentagons and 6 hexagons. Similarly, as 
a result of Stone- Wales type transformation, each C4 de- 
fect region also turns into a nonagon-shaped hole. As 
it was reported experimentally,'^ a hole region is sur- 
rounded by pentagonal and hexagonal rings of C atoms 
to keep the fiatness of the sheet. While the honeycomb 
lattice symmetry of the graphene matrix does not change 
considerably for larger defects Cg and C24, hole region of 
C12 takes almost a circular shape surrounded by regular 
pentagons and hexagons. Apparently, the trend in the 
shape of the hole region is determined by whether the 
edges of domain are zigzag or armchair shaped. 

The energy band structures of graphene nanomeshes 
exhibit interesting variations with respect to their sizes, 
diameters and shapes of nanohole as illustrated in Fig|4] 
(b). While nanomesh of Ci defect have metallic behav- 
ior, periodic structure of C2 defect becomes indirect band 
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FIG. 4: (a) Structural parameters for a nanomesh of C„ hole, 
(b) Large supercells of hexagonal lattice each containing single 
hole of Ci, C2, C4, Ce, C12, C24. The nanomeshes in the third 
row are obtained by saturating Ce and C24 holes by hydrogen 
and also by B and N atoms alternatingly. 



gap semiconductor of 0.65 eV. However, tt and 7r*-bands 
above the Fermi level still very close to each other (0.09 
eV gap) at the iiT-point. C4 defected GNM has metallic 
behavior. The situation becomes even more remarkable 
for larger defects Cg, C12 and C24. GNMs including ei- 
ther Ce or C24 hole have zigzag edges and are metals 
with antiferromagnetic (AFM) ground state. We note 
that the defect induced flat electronic bands around the 
Fermi level occur if the electron spins become unpaired. 
On the other hand, the reconstruction or the dimerization 
of carbon dangling bonds around the defect can cause 
the flat bands in the gap to disappear. Interestingly, 
nanomeshes of Cg and C24 holes, which are metallic, be- 
come semimetal upon the saturation of dangling bonds 
of carbon atoms around the hole by hydrogen atoms. 
This metal-semimetal transformation can have applica- 
tions in graphene-based nanoelectronics. On the other 
hand, GNM of Cg turns semiconducting with a band gap 
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FIG. 5: Band structure of nanomeshes of C12 forming in the 
(nxn) supercells of graphene with n=4...11. 



of 0.1 eV upon the termination of the dangling bonds of 
C atoms around the hole by B and N alternatingly to 
from a B3N3 hexagon. The band opening is explained by 
the breaking of sixfold symmetry due to B3N3 hexagon. 
The formation extended B„N„ honeycomb structure can 
be achieved directly in the course of epitaxial growth of 
graphene and single layer BN.'*^ 

As for GNM with C12, it is a nonmagnetic semimetal 
because of carbon atoms at the edge are dimerized. The 
analysis of the orbital character of linearly crossing tt- 
and 7r*-bands near i^-point using isosurface charge den- 
sities suggests that these bands originate from bonding 
and antibonding combinations of tt- orbitals at the neck 
and around Cg hole. Here we discuss an important as- 
pect of GNMs with C12, that the size of the hexagonal 
supercell or repeat periodicity of C12 is crucial for the 
resulting electronic structure. Fig. [5] shows band struc- 
tures of GNMs including single C12 hole in the supercell 
of (nxn) where n=4...11. For n=4, GNM has the neck 
region consisting of single hexagon is a semimetal. For 
n=5, GNM is a metal and has a neck region which is rel- 
atively ticker, but its C-C bond angles strongly deviate 
from 120°. Surprisingly, GNM with n=6 is a semicon- 
ductor having 1.3 eV band gap. The bond angles still 
continue to deviate from those of graphene. However, 
from n=l on the bond angles at the neck region start 
to be graphene like with regular honeycomb structure. 
Both GNMs with n=7 and 8 are semimetals and have tt- 
and 7r*-bands which linearly cross at Fermi level at the 
if -point. Isosurface charge densities of these bands near 
ii'-points demonstrate that they, in fact, originate from 
the combination of graphene 7r-orbitals. As n increases. 



GNM exhibit the similar trend for n=6-9; namely it is 
semiconductor for n=9, but semimetals for n=10 and 11. 
This variation of band gaps is reminiscent of the family 
behavior of graphene nanoribbons and is related to the 
variation of the thickness of necks between periodically 
repeating C12 holes !2215fl! Here, even if the six fold rota- 
tion symmetry is conserved, the band gap opens for every 
n=3xiV with TV being an integer > 2. However, this gap 
becomes smaller and eventually is closed as n — cx). We 
also note that the family like behavior of GNMs is related 
with the edge structure of the hole. In regard to the size 
of GNM of C12, we note also that the lattice constants 
of corresponding (nxn) supercell is modified with size. 
For example, for n=ll, the lattice constant of supercell 
is contracted by 1%, the contraction is 4% for n=5 and 
40% for n=4. 

Finally, in addition to triangular defect patterns, we 
also discuss the electronic structure of holes arranged in 
a rectangular lattice. In Fig. [6] we show the electronic 
band structures of C12 nanomeshes realized by the su- 
percells of (3x6), (4x8) and (5x10). While C12 holes in 
a small rectangular supercells with small repeat period- 
icity (leading to significant coupling) become semicon- 
ducting, the semimetallic nature indigenous to graphene 
is achieved in large supercells. As shown in Fig. ^c), 
even if the rotation symmetry required for band cross- 
ing is absent, graphene-like Bloch wave functions in the 
rectangular mesh of sparse patterns of C12 holes show a 
semimetallic behavior. This indicates that as the size of 
the supercell becomes larger and the neck gets wider rel- 
ative to the size of the hole, the symmetry requirement 
necessary for the linear band crossing can be relaxed. 



V. MECHANICAL PROPERTIES OF 
NANOMESHES 

Honeycomb structure with sp^ bonding underlies the 
unusual mechanical properties providing very high in- 
plane strength, but transversal flexibility. Here we in- 
vestigate how the mechanical properties of nanomeshes 
generated with patterns of adatoms or holes. We focused 
on the harmonic range of the elastic deformation, where 
the structure responded to strain e linearly. Here e the 
elongation per unit length. The strain energy is defined 
as Eg = Exie) — Et{€ = 0); namely, the total energy at a 
given strain e minus the total energy at zero strain. Nor- 
mally, the Young's modulus is the value, which character- 
izes the mechanical strength of a bulk material. Owing 
to ambiguities in defining the Young's modulus of two 
dimensional structures like a GNM, one can use in-plane 
stiffness C = (I/Aq) ■ (d'^Es/de ^) in terms of the equi- 
librium area of the supercell, AgW^We calculated the 
in-plane stiffness of graphene, and nanomeshes consisting 
of C32H6 , BN substituted graphene (i.e C26[(BN)3] and 
Cq hole in graphene using (4x4) supercell in Fig. [s] The 
stiffness of bare graphene is calculated to be 334 N/m, 
which is in good agreement with the experimental value 
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FIG. 6; (a)-(c) Rectangular patterns of C12 holes repeated in 
(3x6), (4x8) and (5x10) graphene supercells. Atomic struc- 
ture of nanomeshes are given by inset, (d) Semimetallic elec- 
tronic structure and isosurface charge densities of valence (V) 
and conduction (C) bands of (5x10) nanomesh are also shown. 



of 340±50 N/m. Furthermore, the in-plane stiffness val- 
ues of nanomeshes generated on graphene through B3N3 
substitution, Hg adatom pattern and Cg hole are 308 
N/m, 283 N/m and 167 N/m, respectively. Apparently, 
the bare graphene matrix is weakened by the formation of 
any of these nanomeshes. In addition to the calculation 
of in-plane stiffness, we extend our analysis to include 
the plastic deformation region, where the honeycomb like 
structure is destroyed after the yielding point (i.e. on set 
of plastic deformation), and GNM undergoes a massive 
structural deformation. Our preliminary simulations in- 
dicate that the yielding strain of Cg hole GNM is signifi- 
cantly lower than the yielding strains of both C2g[(BN)3] 
and C32Hg GNMs. 



VI. CONCLUSIONS 

While graphene and its various derivatives, GO, CH 
and CF, are important nanomaterials with diverse elec- 
tronic, magnetic and mechanical properties, their proper- 
ties can be modified and multiplied using different meth- 
ods of functionalization. The most pronouncing prop- 
erty of graphene, namely linearly crossing bands at the 
Fermi level and electron-hole symmetry arising there- 
from is usually destroyed, when graphene is functional- 
ized through dopant or vacancy defects. In this work 
we demonstrated that the electron-hole symmetry, even 
Dirac Fermion behavior can be recovered for period- 
ically repeating superstructures or nanomeshes having 
special point group symmetry. In this study we con- 
sidered nanomeshes, which are generated by the deco- 
ration of adatoms, adatom groups or holes, which re- 
peats periodically in graphene matrix. We found that 
types of adatoms and their patterns, the geometry of 
the holes of carbon atoms, the sizes and lattice symme- 
tries of nanomesh provide us with several parameters to 
engineer the electronic and magnetic properties of the 
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nanoniesh. In particular, we showed that by varying 
only the size of the nanomesh including a specific hole 
one can tune between metallic and semiconducting state 
including semimetal with linearly crossing bands. This is 
reminiscent of the family behavior of graphene nanorib- 
bons. 
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